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Abstract. We investigate current fluctuations in non-degenerate semiconductors, on 
length scales intermediate between the elastic and inelastic mean free paths. The shot- 
noise power P is suppressed below the Poisson value Pp isson = 2eJ (at mean current 
7) by the Coulomb repulsion of the carriers. We consider a power-law dependence of 
the elastic scattering time r oc e a on kinetic energy e and present an exact solution of 
the non-linear kinetic equations in the regime of space-charge limited conduction. The 
ratio P/Ppoisson decreases from 0.38 to in the range — | < a < 1. 



1 Introduction 

The noise power P of current fluctuations in an electron gas in thermal equi- 
librium (at temperature T) is related by the Johnson-Nyquist formula P — 
AkTG (with k Boltzmann's constant) to the linear-response conductance G = 
limy-^o dl/dV (with I the mean current in response to an applied voltage V) . 
This formula can be generalized to a large applied voltage, P = 4kT(V/ I)(dl /dV) 2 , 
provided the electron gas remains in local equilibrium with the lattice. Local 
equilibrium requires inelastic scattering. When the conductor is shorter than 
the inelastic mean free path U n and the potential drop V is large enough, the 
Johnson-Nyquist formula no longer applies and a measurement of current noise 
(then also called shot noise) reveals more detailed information about the trans- 
port of charge carriers — in particular about their correlations. The maximal 
noise level Pp isson = 2e/ is attained in absence of all correlations (both in the 
injection process as well as in the subsequent transport). Examples are vacuum 
diodes at large bias in absence of space-charge effects and tunneling diodes with 
low transmissivity. 

Here we consider the transport through a disordered semiconductor of length 
L terminated by two metal contacts, under the conditions of elastic scattering 
(I <C L -C hn, with I the elastic mean free path). In a degenerate conductor 
correlations are induced by the Pauli exclusion principle (for a review of the 
theory of shot noise in this situation see Ref. [1]) and the shot noise has the 
universal value P — |ftoisson [2], [3]. 
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At low carrier concentration the electron gas is non-degenerate, and the Pauli 
principle is ineffective. Because carriers can now accumulate, giving rise to space- 
charge effects, they become correlated through Coulomb repulsion. This is the 
situation which we want to study presently. In a recent Monte-Carlo simulation 
[4] a shot-noise suppression factor of about P/Pp isson = 1/3 was found in the 
regime of space-charge limited transport; an energy- independent elastic scat- 
tering rate was assumed. The coincidence with the noise level obtained in the 
degenerate situation attracted a lot of attention [5]. The degree of universality 
is less pronounced here since the number actually depends on the geometry and 
dimensionality — as well as the scattering mechanism [6] , [7] , [8] . 

In Ref. [6] the problem was investigated for an energy-independent elastic 
scattering time r, using the kinetic theory of non-equilibrium fluctuations (re- 
viewed in Ref. [9]). The non-linear kinetic equations were solved in a certain ap- 
proximation (the drift approximation), with the result P/Pp isson = 0.3410. In 
Ref. [7] we obtained an exact solution, giving -P/-Pp isson = 0.3097, and also con- 
sidered a power-law dependence r ~ e a on the kinetic energy e. For a = — \ (cor- 
responding to short-range impurity scattering or quasi-elastic acoustic phonon 
scattering [10]) we found the exact result P/Pp isson = 0.3777. For other values 
of a we only presented results within the drift approximation. In this work we 
derive the exact solution in the range — | < a < 1. As we will discuss, a should 
be in this range for space-charge limited conduction to be realized. 

2 The drift-diffusion equation 

We consider a three-dimensional conductor of length L and cross-sectional area 
A terminated by two contacts. The equilibrium density p cq of charge carriers 
(charge e, effective mass m) in the decoupled conductor is assumed to be much 
lower than the density p c of those carriers that are energetically allowed (at a 
given voltage V) to enter the conductor from the contacts. (A possible realization 
would be an intrinsic or barely doped semiconductor between two metal contacts 
or two heavily doped semiconducting regions.) The dielectric constant of the 
conductor is k. The temperature T is assumed to be so high that the electron 
gas is degenerate, and a large voltage drop V 3> kT/e is maintained between 
the contacts. Transport is assumed to be diffusive and elastic, I < L < Z; n . We 
assume a power-law energy dependence 



of the elastic scattering time on the kinetic energy e. We want to calculate the 
zero-frequency component 



r(e) = t £' 



.cx 



(1) 




(2) 



of the noise power of the fluctuations SI(t) of the electric current I(t) = I + SI(t) 
around its mean /. 
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We use Cartesian coordinates x, y, z with x parallel to the conductor (the 
current source is at x = 0, the drain at x = L). To linear order in the fluctuations, 
the transverse coordinates can be ignored. In the zero-frequency limit the current 
is independent on x because of the continuity equation and is given by the drift- 
diffusion equation [6], [7] 

I(t) = -J^ J deD(e)p(x,e,t) + E(x,t) J ' de T{x,e,t)^^- + 6J(x,t). (3) 

The electric field E(x,t) is related to the laterally integrated charge density 
p(x, t) by the Poisson equation 

K^E(x,t) = jp(x,t), (4) 

where we omitted the low background charge density — p cq . The fluctuating 
source SJ(x, e, t) accounts for the stochasticity of individual scattering events 
and has the correlator 



6J{x,t)6J(x',t') = 2AS(t-t')6(x-x') J de a(e)J 7 (x,e). (5) 

Here and in Eq. (3), T{x, e, t) = p(x, e, t)/ev(e) with the density of states v(e) = 
4nm(2me) 1 ^ 2 = v$e x l 2 (we set Planck's constant h = 1). The conductivity 
cr(e) = e 2 v(e)D(e) = a e a+3 ^ 2 is the product of the density of states and the 
diffusion constant D(e) = v 2 t/3 = D s a+1 . 



3 Space-charge limited conduction 

For a large voltage drop V between the two metal contacts and a high carrier 
density p c in the contacts, the charge injected into the semiconductor is much 
higher than the equilibrium charge p cq , which can then be neglected. For suf- 
ficiently high V and p c the system enters the regime of space-charge limited 
conduction [11], defined by the boundary condition 

E(x,t)=0 at x = 0. (6) 

Eq. (6) states that the space charge Q = p(x) dx in the semiconductor is 
precisely balanced by the surface charge at the current drain. At the drain we 
have the absorbing boundary condition 

p(x,t) =0 at x = L. (7) 

With this boundary condition we again neglect p oq . 

To determine the electric field inside the semiconductor we proceed as follows. 
Since scattering is elastic, the total energy u = e — e<f>(x, t) of each carrier is 
preserved. The potential gain —e<f>(x, t) (with E = —d<t>/dx) dominates over 
the initial thermal excitation energy of order kT almost throughout the whole 
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semiconductor; only close to the current source (in a thin boundary layer) this 
is not the case. We can therefore approximate the kinetic energy e w — e<j> and 
introduce this into D(s) and da/de. Substituting into Eq. (3) one obtains 



<7 [-e<f){x' ,t)] a + 3 / 2 '' 



D J x [-e4>{x',t)r+^ 

(9) 

where the absorbing boundary conditions have been used. From the Poisson 
equation (4) we find the third-order, non-linear, inhomogcneous differential equa- 
tion 

2(-<£)<W' + 4(^) Q +V" = BI[1 + 6i(x, t)], (10) 
6i{ x ,t)= I ® (11) 



for the potential profile <j>(x, t). Primes denote differentiation with respect to x, 
and B = 6/e a /j,onA with /io = ero/m. 

Since the potential difference V between source and drain does not fluctuate, 
we have the two boundary conditions 0(0, t) — 0, <j>(L, t) = —V. Eqs. (6) and (7) 
imply two additional boundary conditions, (p'(0,t) — 0, <f>"(L,t) = 0. 

The differential equation (10) and the accompanying boundary conditions 
possess two remarkable scaling properties: The product BI of material parame- 
ters and mean current I and the length L can be eliminated by introduction of 
the scaled potential 

X (x, t) = - (L 3 BI)- 1/(a+2) <j>{xL, t). (12) 
The rescaled differential equation reads 

2xVx" - 4 X Q + V" = l + 6i, (13) 

which has to be solved with the boundary conditions x(0,t) = 0, x(l,£) = 

(L 3 Biy 1/{a+2) V, x'(0,i) = 0, x"(M) = 0. The scaling properties entail that 
the shot-noise suppression factor depends only on the exponent a, but no longer 
on the parameters L, A, V, t , and n. 

We will solve this boundary value problem for x — X + ^Xi first f° r t h e rnean 
(Section 4) and then for the fluctuations (Section 5), in both cases neglecting 
terms quadratic in 5\- 
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4 Average profiles and the current-voltage characteristic 

The averaged equation (13) for the rescaled mean potential x(x) reads 



We seek a solution which fulfills the three boundary conditions x(0) = 0, x'(0) = 
0, x"(l) = 0. The value of x at the current drain determines the current-voltage 
characteristic 



We now construct x( x )- The function Xo( x ) = ao^ with (3 = 3/(2 + a) and 
a = [2/3(/3 — 1)(4 — /3)]~^ 3 solves the differential equation and satisfies the 
boundary conditions at x = 0, but x'o( x ) for any finite x. We substitute into 
Eq. (14) the ansatz x(x) — a i xll+ ^ \ consisting of Xo( x ) times a power series 
in x 7 , with 7 a positive power to be determined. This ansatz proves fruitful since 
both terms on the left-hand side of Eq. (14) give the same powers of x, starting 
with order x° in coincidence with the right-hand side. By power matching one 
obtains in first order the value for ao given above. The second order leaves a\ as 
a free coefficient, but fixes the power 7 = (8 - 5/3 + ^-32 + 40/3 + /3 2 )/4 . The 
coefficients a; for I > 2 are then given recursively as a function of a\, which is 
finally determined from the condition x"(l) = 0. 

In Fig. 1 the profiles of the potential <f> oc x\ the electric field E oc X'i an d 
the charge density p oc x" ar e plotted for various values of a. The coefficient 
X(l) appearing in the current-voltage characteristic (15) can be read off from 
this plot. The behavior at the current source changes qualitatively at a = — -| 
(see Section 7). 

5 Fluctuations 

The rescaled fluctuations Sx{x, t) = "4>(x, t) fulfill the linear differential equation 



The solution of the inhomogeneous equation is found with help of the three 
independent solutions of the homogeneous equation £[ip] = 0, ipi{x) — x'( x ), 
ip2(x) = x(x) - (x/(3)x'(x), and 



2X Q X'X"-4X Q+1 X"' = 1. 



(14) 




(15) 



C[il>\ = -A X a+ V" + 2x a xV" + 2x"x"V / 

+ [2ax Q ~ 1 x'x" - 4(a + l)x Q x"'] V> = 



(16) 




(17) 
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x/L 



Fig. 1. Profile of the mean electrical potential <f> [in units of (L 3 BI) 1 '( a+2 ' , with 
B = 6m/e a+1 T nA], the electric field E [in units of (L A BI) 1/(a+2) / L], and the charge 
density p [in units of k(L 3 B/) 1 /( q: + 2 ) /L 2 ], following from Eq. (14) for different values 
of a. 
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where we have defined W{x) = ^pi{x)^p' 2 {x) — ip' 1 (x)tp2(x) . The special solution 
which fulfills V(0, t) = ip'(0, t) = t) = is 



V>(x, t) = J dx' [0(x - x')Vi(x)V 2 (x') + 0(x' - x)M*)Mx) 



^i(l) (| , w ,,ir, „ 6I(t)-6J(x»,t) W(x") 



The condition ip"(l,t) — relates the fluctuating current 5/ to the Langevin 
current S J. The resulting expression is of the form 

5I(t)=C- 1 [ dx5J(x,t)g(x), (19) 
Jo 

with the definitions C = Jq dx £ (x) , 

W(x) / (1 - l//3)x /2 (l) f 1 x 1/2 ^)M^) \ 

f w - ^ + 3 /2(x) ^ + 4^+1/2(1)^(1) j x w 2 (xo ; • (2Uj 

The shot-noise power is found by substituting Eq. (19) into Eq. (2) and in- 
voking the correlator (5) for the Langevin current. This results in 

P = 2 ( L dx (^) H(x) (21) 



C 

with H(x) = 2AJ dea(e)T(x,e) w 2cr o [-e0(x)] a+3 / 2 ^(x). Eq. (8) gives 

H(x) = 2e/x Q+3/2 d) £ dx' - a+3 ) 2{xl) = 4P Poisson x a+1 d)x"d), (22) 

where we integrated with help of Eq. (14) and used x"(l) = 0. 

In Fig. 2 we plot the ratio P/ Pp isson as a function of the parameter a (solid 
curve). The shot-noise suppression factor P/Ppoisson = 0.3777 for a = — | and 
goes to zero as a — > 1. 



6 Drift approximation 

A simple formula for the shot-noise suppression factor can be found when one 
neglects the diffusion term in Eq. (3) and considers instead of Eq. (13) the corre- 
sponding differential equation (4a + 6)x Q x'x" = 1 + Si . The is the drift approx- 
imation of Ref. [6]. The order of the differential equation is reduced by one, so 
that we also have to drop one of the boundary conditions. The absorbing bound- 
ary condition x"(l, t) = is the most reasonable candidate, because even for the 
resulting mean profile \{x) = b^x 13 with (3 = 3(2 + a) _1 and b = [f3 2 (j3- l)]^/ 3 
most carriers remain concentrated close to the current source. The differential 
equation for the fluctuations aip/x + tp'/x' + i ) " /x" — & can be solved with 
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-0.5 -0.25 0.25 0.5 0.75 1 
a 

Fig. 2. Shot-noise power P as a function of a. The exact result (solid curve) is compared 
with the approximate result (24) (dashed curve) . 



help of the homogeneous solutions ?pi(x) — x@ 1 and ip2(x) — x 3 2 @ . The inho- 
mogeneous solution that fulfills -0(0, t) = 0, tp'(0, t) = is 

^ t ) = ~p f da;' [x 3 - 2 ^' 3 ' 3 - 4 - x^ 1 ] <5i(z', t). (23) 

We demand that the voltage does not fluctuate, ipO-y t) = 0, and obtain Eq. (19) 
with now G(x) = 1 — x 3 P~ i . The shot noise power is finally found from Eq. (21) 
with H{x) = P Poisson x 3 ^/ 2 Jl dx> x<Pli-\ 

_ 6(a-l)( a + 2)(16a 2 + 36a~157) 
7 Poisson ~ 5(2a-5)(8a-17)(13 + 8a) ' 1 j 

This is the dashed curve in Fig. 2. 



7 Discussion 

The shot-noise suppression factor P/Pp isson varies from 0.38 to in the range 
— \ < a < 1, which includes the case of an energy- independent elastic scattering 
rate (a = 0, P/Pp isson = 0.3097) and the case of short-range scattering by 
uncharged impurities or quasi-elastic scattering by acoustic phonons (a = — ^, 
P/Ppoisson = 0.3777). The results in the drift approximation (24) are about 
10% larger. Our values are somewhat smaller than those following from the 
numerical simulations of Gonzalez et al, who found P/Pp isson = | for a = 
[4] and P/Ppoisson - 0.42 - 0.44 for a = -\ [12]. 
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Our considerations require the exponent a to be in the range — ^ < a < 1. 
For a < — \ the mean free path I oc e a+1 / 2 diverges at small kinetic energies. 
The carriers at the current source therefore enter the conductor ballistically and 
accumulate only at a finite distance from the injection point. Fig. 1 indicates 
that the charge density at the current source must be zero if one insists that the 
electric field vanishes. Nagaev [8] has shown that full shot noise, P = -Ppoisson, 
follows for a = — §. Presumably, -P/-Pp isson will decrease monotonically from 1 
for a = — | to 0.38 for a = — |, but wc have no theory for this range of en's. 
For a > 1 the resistance R becomes infinitely large, because the coefficient x(l) 
in the current-voltage characteristic (15) diverges. An intuitive understanding 
can be obtained by equating the potential gain <j> ~ (Dt) 3 ^ 2a+4 ^ (acquired by 
diffusing close to the current source for a time t) with the increase in kinetic 
energy e: For a > 1 this time t oc e( 1_Q! )/ 3 is seen to diverge for small e. We 
found that the shot-noise power vanishes as a — ► 1. Presumably, a non-zero 
answer for P would follow for a > 1 if the non-zero thermal energy and finite 
charge density at the current source is accounted for. This remains an open 
problem. 

Discussions with O. M. Bulashenko, T. Gonzalez, J. M. J. van Leeuwen, and 
W. van Saarloos are gratefully acknowledged. This work was supported by the 
European Community (Program for the Training and Mobility of Researchers) 
and by the Dutch Science Foundation NWO/FOM. 
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